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Abstract. We consider iV-wave type equations related to the orthogonal algebras obtained 
from the generic ones via additional reductions. The first Z2-reduction is the canonical one. 
We impose a second ^-reduction and consider also the combined action of both reductions. 
For all three types of iV-wave equations we construct the soliton solutions by appropriately 
modifying the Zakharov-Shabat dressing method. We also briefly discuss the different types 
of one-soliton solutions. Especially rich are the types of one-soliton solutions in the case when 
£h ' both reductions are applied. This is due to the fact that we have two different configurations 

of eigenvalues for the Lax operator L: doublets, which consist of pairs of purely imaginary 
eigenvalues, and quadruplets. Such situation is analogous to the one encountered in the 
sine-Gordon case, which allows two types of solitons: kinks and breathers. A new physical 
system, describing Stokes-anti Stokes Raman scattering is obtained. It is represented by 
a 4-wave equation related to the B2 algebra with a canonical Z2 reduction. 
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1 Introduction 

The iV-wave equation related to a semisimple Lie algebra q is a matrix system of nonlinear 
differential equations of the type 



H| i[J,Q t (x,t)]-i[I,Q x (x,t)] + [[I,Q(x,t)],[J,Q(x,t)]]=0, (1.1) 

where the squared brackets denote the commutator of matrices and the subscript means a partial 
derivative with respect to independent variables x and t. The constant matrices I and J are 
regular and real elements of the Cartan subalgebra b, of the Lie algebra q. The matrix- valued 
function Q(x, t) G fl can be expanded as follows 

Q(x,t) = ^2 Q a (x,t)E ai 

where A denotes the root system and E a are elements of the Weyl basis of q parametrized 
by roots of q. It is also assumed that Q(x,t) satisfies a vanishing boundary condition, i.e. 
lim Q(x,t) = 0. 

x— >±oo 
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The iV-wave equation is an example of an S-integrable evolutionary equation. Such type 
of equations appear in nonlinear optics and describes the propagation of N wave packets in 
nonlinear media [2]. 

Another application of the N-w&ve systems is in differential geometry. Ferapontov [3] showed 
that N-w&ve equations naturally occurred when one studied isoparametric hypersurfaces in 
spheres. 

Our aim in this paper is two-fold. First, we outline the derivation of N-w&ve equations with 
Z2-reductions and calculate their soliton solutions by applying one of the basic methods in theory 
of integrable systems - Zakharov-Shabat dressing procedure. One of the considered examples has 
a physical interpretation - it describes a Stokes-anti-Stokes Raman scattering in nonlinear optics. 
Secondly, we further reduce the N-w&ve equations by imposing a second ^-reduction. As a result 
we derive a special class of N-w&ve equations which, like the sine-Gordon equation (SG), possess 
breather type solutions. Next we obtain their soliton solutions. The additional symmetries of 
the nonlinear equations have been taken into account when we chose a proper dressing factor. 
The soliton solutions of a Z2 x ^-reduced N-w&ve equation are of two types. The solitons of 
the first type are connected with a couple of purely imaginary discrete eigenvalues of the Lax 
operator L. These solutions correspond to the kinks (or topological solitons) of SG. The solitons 
of the second type, associated with quadruplets of discrete eigenvalues L which are symmetrically 
located with respect to the coordinate axis in the complex A-plane, are analogues of the breather 
solutions of SG. Recently the Darboux-dressing transformations have been applied to the Lax 
pair associated with systems of coupled nonlinear wave equations: vector nonlinear Schrodinger- 
type equations [3j. Solutions with boomeronic and trapponic behaviour were investigated. 

The problem for classification and investigation of all admissible reductions of an integrable 
equation is one of the fundamental problems in the theory of integrable systems. N-w&ve 
equations with canonical Z2 symmetries have been discussed for the first time by Zakharov and 
Manakov in [23] for q ~ sl(n). More recently the Z2-reductions of the N-w&ve equations related 
to the low-rank simple Lie algebras were analyzed and classified [8j [9] . 

Since we shall deal with the inverse scattering transform we begin with a reminder of all 
necessary facts concerning that theory. For more detailed information we recommend the mono- 
graphs [21] and [21]. 

2 General formalism 

As we mentioned above the N-w&ve system (jl.ip is an integrable one. It admits a Lax repre- 
sentation with Lax operators 

Ltp(x, t, A) = (id x + U(x, t, A)) ip(x, t, A) = 0, 



where A is an auxiliary (so-called spectral) parameter and the potentials U(x,t,X), V(x,t, A) 
are elements of and they are linear functions of A defined by 



The nonlinear evolution equation itself is equivalent to the compatibility condition of the diffe- 
rential operators L and M 



Mij}(x, t, A) = (id t + V(x, t, A)) ip(x, t, A) 




U(x, t, A) = U°(x, t)-\J = [J, Q(x, t)] - AJ, 
V(x, t, A) = V°{x, t)-XI = [I, Q(x, t)] - XI. 



[L,M]=0 o i[J,Q(x,t) t ]-i[I,Q(x,t) x ] + [[I,Q(x,t)],[J,Q(x,t)]]=0. 
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Since L and M commute they have the same eigenfunctions and the system (|2.ip can be pre- 
sented by 

Ltp(x, t, A) = (id x + U(x, t, A)) ip(x, t, A) = 0, 

Mip(x, t, A) = (id t + V(x, t, A)) if>(x, t, A) = if>(x, t, A)C(A), (2.2) 

where C(A) is a constant matrix with respect to x and t. The fundamental solutions tp(x,t, A) 
of the auxiliary linear system f|2.2 j) take values in the Lie group G which corresponds to the Lie 
algebra q. 

In order to find the soliton solutions we need the so-called fundamental analytic solutions 
(FAS) derived for the first time in [23] for G ~ SL(n); for other simple Lie algebras see [7j. 
There is a standard algorithm to construct these solutions by using another class of fundamen- 
tal solutions of the linear problem (I2.2D - Jost solutions. The Jost solutions ijj±(x,t, A) are 
determined by their asymptotics at infinity, i.e. 

lim ii±(x,t,\)e iXJx = 1. 
£^±00 

Remark. This definition is correct provided we fixed up the matrix- valued function C(A) by 
C(A) = lim V(x,t,X) = -XI, 

x— ► ±oo 

i.e. the asymptotics of ip±(x, t, A) are i-independent. C(A) is directly related to the dispersion 
law of the nonlinear equation. Thus the dispersion law of the iV-wave equation is a linear 
function of the spectral parameter A. 

The Jost solutions ip±(x,t,X) are linearly related, which means that there exists a matrix 
T(t, A) such that 

i/>-(x,t,\) = il>+(x,t,\)T(t,\). 

The matrix-valued function T(t, A) is called a scattering matrix. Its time evolution is determined 
by the second equation of (12.2|) . i.e. 

id t T(t,X)-\[I,T(t 1 X)]=0. 

Consequently 

T(t, A) = e" iA/ *r(0, X)e iXIt . (2.3) 

The inverse scattering transform (1ST) allows one to solve the Cauchy problem for the nonli- 
near evolution equation, i.e. finding a solution Q(x, t) when its initial condition Q- m {x) = Q(x, 0) 
is given. The idea of 1ST is illustrated in the following diagram 

Q in (x) -► U(x, 0, A) ^ T(0, A) — ► T(t, A) ^ U(x, t, A) -» Q(x, t). 

The first step consists in constructing the scattering matrix at some initial moment t = 
by using the potential at the same moment (or equivalently by the solution Qi n (x) at that 
moment). This is a direct scattering problem (DSP). The evolution of the scattering matrix 
(data) is already known and it is given by (|2.3p . The third step is recovering of the potential 
U (x, t, A) and respectively the solution of the nonlinear equation Q(x, t) at an arbitrary moment 
from the scattering data at that moment - this is an inverse scattering problem (ISM). That step 
is actually the only nontrivial one. Thus following all steps we can solve the Cauchy problem for 
the nonlinear evolution equation. Since we know the evolution of scattering data we can easily 
determine the time dependence of fundamental solutions, solutions of nonlinear problem etc. 
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The Jost solutions are defined for real values of A only (they do not possess analytic properties 
for A ^ R). For our purpose it is necessary to construct fundamental solutions which admit 
analytic continuation beyond the real axes. It can be shown that there exist fundamental 
solutions x + ( x ,t,X) and x~( x ,t,X) analytic in the upper half-plane C + and in the lower half- 
plane C_ of the spectral parameter respectively. They can be obtained from the Jost solutions 
by a simple algebraic procedure proposed by Shabat [TU [20] , see also [24] . The procedure uses 
a Gauss decomposition of the scattering matrix T(t, A), namely 

X ± (x, t, A) = V_(x, t, X)S ± (t, A) = t, X)T*(t, X)D±(X), 

where matrices 5 (t, A), T (t, A) and D ± (X) are Gauss factors of the matrix T(t, A), i.e. 

T(t,\) = TT(t,\)D ± (\)(S ± (t,\)r 1 . 

The matrices S + (t, A) and T + (t,X) (resp. S~(t, A) and T~(t,X)) are upper (resp. lower) trian- 
gular with unit diagonal elements and taking values in G. Their time dependence is given by: 

id t S±(t, A) - X[I, S ± (t, A)] = 0, idt^it, A) - A[I, r±(t, A)] = 0. 

The matrices D + (X) and D _ (A) are diagonal and allow analytic extension in A for ImA > and 
ImA < 0. They do not depend on time and provide the generating functionals of the integrals 
of motion of the nonlinear evolution equation pH EJ E], [9] , see also the review paper [7] . 

A powerful method for obtaining solutions to nonlinear differential equations is Backlund 
transformation (see [15] and [181117] for more detailed information). A Backlund transformation 
maps a solution of a differential equation into a solution of another differential equation. If both 
equations coincide one speaks of an auto-Backlund transformation. A very important particular 
case of an auto-Backlund transformation is the dressing method proposed by Zakharov and 
Shabat [25]. Its basic idea consists in constructing a new solution Q(x, t) starting from a known 
solution Qo(x,t) taking into account the existence of the auxiliary linear system (|2.2p . 

Let ipo(x,t,\) satisfy the linear problem 

L ^ {x, t, A) = id x M%, t, ^ + ([J, Qo(x, t)] - \J)M%, t, A) = 0. (2.4) 

We construct a function vjj(x,t,\) by introducing a gauge transformation g(x,t,X) - dressing 
procedure of the solution ipo(x, t, A) 

rpo(x, t, A) -> ip(x, t, A) = g(x, t, \)ipo{x, t, A) 



such that the linear system (12. 4p is covariant under the action of that gauge transformation. 
Thus the dressing factor has to satisfy 

id x g(x,t,\) + [J,Q(x,t)]g(x,t,\) - g(x,t,X)[J,Q (x,t)] - X[J,g(x,t, X)] =0. (2.5) 

If we choose a dressing factor which is a meromorphic function of the spectral parameter A as 
follows 

A(x,t) B(x,t) , s 

g(x,t,X) = Z + ^-j± + ^-jl, (2.6) 

where A ± G C±, we obtain the following relation between Q(x,t) and Qo(x,t) 
[J,Q(x,t)} = [J,Q (x,t) + A(x,t) + B(x,t)}. 



As a result we are able to construct new solutions if we know the functions A(x,t) and B(x,t). 
We will show later how this can be done. 
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The simplest class of solutions are the so-called reflectionless potentials. A soliton solution 
is obtained by dressing the trivial solution Qo(x,t) = 0. Then the fundamental solution of the 
linear problem is just a plane wave ip (x,t,X) = e - iX( - Jx+It 1 and the one-soliton solution itself is 
given by 

[J, Qn (x, t)] = [J, A u (x,t) + B u (x,t)}. 

As we said above the dressing procedure maps a solution of the linear problem to another so- 
lution of a linear problem with a different potential. In particular the Jost solutions ipo±(x, t, A) 
are transformed into 

ip±(x,t,\) = g(x,t,\)Tp 0t ±{x,t,\)g± 1 (\), 

where the factors g± (A) ensure the proper asymptotics of the dressed solutions and are defined by 

g±(X) = lim g(x,t, A). 

x— >±oo 

Hence the dressed scattering matrix T(t, A) reads 

T(t,X)=g + (X)T (t,X)gZ 1 (X). 
It can be proven that the FAS Xo( x ,t, A) transform into 

X ± (x,t,X) = g(x,t,X) X ^(x,t,X)gZ 1 (X). (2.7) 

The spectral properties of L are determined by the behaviour of its resolvent operator. The 
resolvent R(t, A) is an integral operator (see [U [7] for more details) given by 



R(t,X)f(x) = / TZ(x,y,t,X)f(y)dy, 
J —00 

where the kernel 1Z(x, y, t, A) must be a piece-wise analytic function of A satisfying the equation 

LTZ{x, y, t, A) = S(x — y)1. 

The kernel lZ(x, y, t, A) can be constructed by using the fundamental analytical solutions as 
follows 



K(x,y,t,X) 



TZ + (x,y,t, A), ImA > 0, 
1Z~(x, y, t, A), ImA < 0, 



where 



^(x, y, t, A) = ±i X ± (x, t, X)B ± (x - y)(x ± (y, A))" 1 , 

e±( x - y ) = e( T (x - y))u - 0(±(x - y ))(i - n), 

II — ^ ] Epp , {,Epq) rs — 5p r 5q S , 

p=l 

where 9 is the standard Heaviside function and a is the number of positive eigenvalues of J. 
Due to the fact that we have chosen J to be a real matrix with 

Jl > J 2 > • • • > Ja > > J a +1 > ■■■> J n , 

the resolvent R(t, A) is a bounded integral operator for ImA 7^ 0. For ImA = R(t, A) is an 
unbounded integral operator, which means that the continuous spectrum of L fills up the real 
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axes M. Since the discrete part of the spectrum of L is determined by the poles of R(t, A) it 
coincides with the poles and zeroes of x^{ x ^^)- 

From (|2.7|) and from the explicit form of 1Z(t, A) it follows that the dressed kernel is related 
to the bare one by 

^(z.y,*, A) = g(x,t,\)KQ(x,y,t,\)g~ 1 (y,t,\). 

If we assume that the "bare" operator Lq has no discrete eigenvalues then the poles of g(x, t, A) 
determine the discrete eigenvalues of L 

Lq — — > L 44> spec(Lo) — ► spec(L) = spec(Lo) U {A + , A - } . 

Many classical integrable systems correspond to Lax operators with potentials possessing 
additional symmetries. That is why it is of particular interest to consider the case when certain 
symmetries are imposed on the potential U(x,t, A) (resp. on the solution Q(x,t)). 

Let Gr be a discrete group acting in G by group automorphisms and in C by conformal 
transformations 

k : A — > k(A). 

Therefore we have an induced action of Gr in the space of functions f(x, t, A) taking values in G 
as follows 

JC : f(x,t,X) - f(x,t,X) =K{f (z.t.K-^A))) , K E Aut(G). 

This action in turn induces another action on the differential operators L and M 

L(X) = /CL^^A))/^ 1 , 
M(A) =KM(k~ 1 (\))1C- 1 . 

The group action is consistent with Lax representation which is equivalent to invariance of Lax 
representation under the action of Gr, i.e. if [L, M] = then also 

[L,M] =0. 

The requirement of G/j-invariance of the set of fundamental solutions {ip(x,t, A)} leads to 
the following symmetry condition 

kuix^n^iXfik- 1 = u{ x ,t,\). 

In other words the potential U(x,t,X), as well as Q(x,t) are reduced. This fact motivates the 
name of the group Gr - a reduction group, proposed by Mikhailov in |16j . 

One can prove that the dressing factor g(x, t, A) ought to be invariant under the action of the 
reduction group Gr 

k(g(x,t,K- 1 (X)))=g(x,t,X). (2.8) 
3 AT-wave equations with a Z, 2 reduction 

In this chapter we are going to demonstrate an algorithm to derive soliton solutions for a iV-wave 
equation related to the orthogonal algebra (i.e. q = 5o(n, C) ) with a Z2 reduction imposed on it. 
We pay special attention to the particular case of B2 — so (5, C) algebra: we display the explicit 
form of the 4-wave system and calculate by components its one-soliton solutions. One of the 
4-wave systems under consideration enjoys a physical application in nonlinear optics. We are 
following the ideas presented by Zakharov and Mikhailov in [22]. 
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3.1 Physical applications of 4- wave equation with a Z2 reduction 

From now on we shall deal with the orthogonal group SO(n,C) and its Lie algebra so(n, C). 
That is why we remind the reader several basic facts concerning this topic without proving them, 
for more detailed information we recommend the book |13j . 

The orthogonal group is, by definition, a matrix group which consists of all isometries acting 
on C n , i.e. 

SO{n, C) = {C e GL(n, C)| C T SC = S) , 
where S stands for the metric in C™. It determines a scalar product by the formula 

(u, v) = u T Sv, u, v £ C n . 
so(n, C) is the matrix Lie algebra of all infinitesimal isometries on C n , i.e. 

so(n, C) = {c G g[(n, C)| c T 5 + 5c = 0} . 

The orthogonal algebra is a semisimple Lie algebra. We remind that its maximal commutative 
subalgebra is called a Cartan subalgebra f). It is more convenient to work with a basis in C n 
such that the matrix of S in that basis reads 
n 

S = Y,(-l) k ~ 1 (Ek,2n+i-k + E 2n+l ^ k ) + {-l) n E nn , for B n ~50(2n + l,C) 
fc=i 

and 

n 

S = ^2(-l) k ~ 1 (E k , 2n+l _ k + E 2n+l _ ktk ), for D n ~so(2n,C). 
k=i 

This choice of a basis ensures that the corresponding Cartan subalgebras consist of diagonal 
matrices. More precisely, Cartan basis {H k }^ =1 , i.e. the basis which spans f), has the form 

H k = E kk — E 2n+ 2~k,2n+2-k, f° r B n , 

H k = E kk — i?2n+l-fc,2n+l-fci for D n . 

The root systems of the series B n and D n consist of the roots 

±ei±ej, ei, i^j, for B n , 
±ei±ej, i^j, for D n , 

where {ej}™ =1 stands for the standard orthonormal basis in C n . The positive roots are as follows 

ei±ej, i<j, ei, for B n , 
&i±ej, i < j, for D n . 

The set of all simple roots reads 

Oi k = &k - ejfc+i, k = l,...,n-l, a n = e n for B n , 
a k = e k - e k+ i, a n = e n _i + e n , for D n . 

Let us illustrate these general results by a physical example related to the B2 algebra. This 
algebra has two simple roots a± = e± — e 2 , a 2 = e 2 , and two more positive roots: a\ + a 2 = e± 
and a\ + 2a 2 = e± + e 2 = a max . When they come as indices, e.g. in Q a (x,t) we will replace 
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them by sequences of two positive integers: a — > kn if a = ka\ + na 2 ; if a = — (ka\ + 7102) we 
will use kn. The reduction which extracts the real forms of B2 is 

K x U\x, t, \*)K{ 1 = U{x, t, A), 

where K\ is an element of the Cartan subgroup: K\ = diag (si, s 2 , 1, s 2 , si) with = ±1, 
k = 1,2. This means that Jj = J*, i = 1,2 and Q a (x,t) must satisfy: 

QtoO^ = -si^QwO; Qtt( x ' *) = SiQu(x, t), 

Qj2-(x, t) = -sis 2 Q* 12 { x i *)> Qor^' *) = - s 2QSi( x ' *)• 

and the matrix Q(x,t) is defined by 

= Qio(x,t)E ei _ e2 + Q 12 (x,t)E ei+e2 + Qn(a;,f)£ ei + Qoifc, i)#e 2 
+ Qio(s,t)£_ 

(ei-e 2 ) + Ql2( X ' ^)-^-(ei+e 2 ) + Q\\{ x i t)E- ei + Q 01 (a;, t)E^ e2 . 

As a result we get the following 4-wave system 

i{Ji - J 2 )Qio,t( x ,t) - i(h - h)Qw,x{x,t) - ks 2 Qu(x,t)QQ 1 (x,t) = 0, 

iJ\Qn,t{x,t) - ihQu >x (x,t) - k(Q w Q i + s 2 Q 12 Q* Ql )(x, t) = 0, 
i(Ji + J 2 )Qi 2)t (x,t) - i(h + h)Qi2,x(x,t) - kQ u (x,t)Q 01 (x,t) = 0, 
iJ2Qoi,t( x ' *) ~ ihQoi,x(x, t) - ksi(Q* u Qi 2 + s 2 Q* w Qn)(x, t) = 0, 
where k = J\I 2 — J 2 I\ and its Hamiltonian H(t) = Ho(t) + Hi nt (t) is 

H (t) = - dx [(h - I 2 )(Q 10 (x,t)Q* Wx (x,t) - Q 10tX (x,t)Q* w {x,t)) 
^ J -00 

+ h{Qoi(x, t)Qo l x (x, t) - Qot,x(x, t)Qoi( x > 0) + h(Qu{x, t)Qn, x (x, t) 

- *)<2u(^ *)) + (A + h)(Qi2{x, t)Q* 12jX (x, t) - Q 12 , x (x, t)Q* l2 (x, t))] , 

/oo 
dx [s 2 (Qi 2 (x,t)Q* u (x,t)QQ l (x,t) + Q]; 2 (x,t)(5ii(2;,t)(5oi(2;,*)) 
-00 

+ (Qn(x, t)Q* 01 (x, t)Q* 10 (x, t) + Q* n (x, t)Q 0l (x, t)Q l0 (x, t))} , 
and the symplectic form: 

/oo 
dx [(Ji - J 2 MQio(x,0 A<5<%(x,i) + J 2 8Q i(x,t) A6Q^(x,t) 
-00 

+ Ji(5Qn(x, i) A t) + (Ji + J 2 )-5Qi 2 (x, i) A £Q* 2 (a;, t)] . 

The invariance condition f|2.8[) leads to the following form of the dressing factor 



= 1 + x^aT + a^TF? ' (3 ' 1} 



i.e. comparing with (|2.6j) we see that 

B(x,t) =lfi5A*(x,t)(iifi5)" 1 , A" = (A H 



By taking the limit A — > 00 in equation (|2.5p and taking into account the explicit formula (|3.ip 
one can derive the following relation between the bare solution Qo(x,t) and the dressed one 

Q(x,t) 

[J,Q(x,t)\ = [J,Q (x,t) + A(x,t) + K 1 SA*(x,t)SK 1 ], 
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Thus the one soliton solution of the iV-wave equation is determined by one matrix- valued func- 
tion A(x, t). We will obtain A(x,t) in two steps by deriving certain algebraic and differential 
relations. First of all recall that the dressing factor g(x, t, A) must belong to the orthogonal 
group SO(n, C), hence 

g- 1 (x,t,X) = S- 1 g T (x,t,X)S. 

Besides, the equality gg~ l = 1 must hold identically with respect to A, therefore A(x, t) satisfies 
the algebraic restrictions 

A(x, t)SA T (x, t) = 0, A(x, t)Suj T (x, t) + co(x, t)SA T (x, t) = 0, (3.2) 

where 

, . „ K 1 SA*(x,t)SK 1 . 

u(x,t) = t-\ -, X + = u + iu. 

y ' ' 2*1/ ' 

From the first equality in (|3,2p it follows that the matrix A(x, t) is a degenerate one and it can 
be decomposed 

A(x,t) = X(x,t)F T (x,t), 

where X(x, t) and F(x, t) are n x k (1 < k < n) matrices of maximal rank k. The equalities (|3.2p 
can be rewritten in terms of X(x,t) and F(x,t) as follows 

F T (x, t)SF{x, t) = 0, X(x, t)F T {x, t)Suj T {x, t) + u(x, i)SFX T {x, t) = 

or introducing a k x k skew symmetric matrix a(x,t) the latter restriction reads 

/ + KtSX-^FH^SKA _ 
\ 2iu J 

Another type of restrictions concerning the matrix-valued functions F(x, t) and a(x, t) comes 
form the A-independence of the potential [J, Q(x,t)\. If we express the potential in the equa- 
tion (12. 5|) we get 

[J, Q(x, t)] = -id x gg~ x (x, t, A) + g[J, Q (x, t)} g" 1 (x, t, A) + A ( J - gJg^ 1 (x, t, A)) . (3.4) 
Annihilation of residues in (|3.4p leads to the following linear differential equations 
id x F T {x,t) - F T (x,t)([J,Q (x 1 t)\ -A+J) = 

and 

id x a(x, t) + F T (x, t)JSF(x, t) = 0. 
After integration we obtain that 

F(x,t) =Sxt(x,t,X+)SF , 
where the constant matrix Fq obeys the following equality 

f£sf = 

and the matrix a(x,t) reads 

a(x, t) = FT( x +(x, t, A+))- 1 a AX ^(x, t, X + )SF + a . 
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In the soliton case the fundamental solution is just a plane wave. Therefore functions F(x, t) 
and a(x,t) get the form 

F(x, t) = e iX+{Jx+It) F , a{x, t) = iF^(Jx + It)SF + a . 

It remains to find the factor X(x,t). For this purpose we solve the linear equation (|3.3|) and 
find that 

X(x,t) = 2v(iK 1 F* - 2vSF{F^K 1 F)~ 1 a*)(F T K 1 F* - Av 2 a{F^ RxF)' 1 a*)^ 1 . 

In particular, in the simplest case when rankX(a;,t) = rankF(x,i) = 1 we have a(x,t) = and 
therefore 



X(x,t) 



2iv 



F j !(x,t)K 1 F{x,t) 
The soliton solution obtains the form 

2iv 



KiF*(x,t). 



Qij{x,t) = < 



[(KiF*(x, t))iFj(x, t) 



F^(x,t)K 1 F(x,t) 

+ (-l)^+i J F 6 _ i (x,t))(^ 1 F*(x,t)) 6 _ j ), 
0, i = j. 

Turning back to the B 2 case this result can be written by components as follows 
2%v 



Qw(x,t) 



F^KiF 
2iv 



j 77* 77 Z2 — \ + * zi) i „ 77* 77 j(\ + *Z2 — \ + zi) 



2%v 



-i A + zi+X 



" 22 )+ S2 F * 2 F , 5 e-< A+ ^ +A+%2 )), 
iX+t 



Qi2(x,t) = pYiCF ( SlF o* lFo ' 4e " 

Qoi(x,t) = ( S2J P * 2 Fo, 3 e- a+ ^ 2 + F^F 0A e-^ z > 

F^K X F = Sl \F 0A \ 2 e- 2 ^ + s 2 \F , 2 \ 2 e- 2 ^ + |F , 3 | 3 + s 2 \F 0A \ 2 e 2 ^ + S i|F , 5 |V^\ 
z a = J a x + I a t, a =1,2. 

A natural choice for the matrix which determines the action of the reduction group Gr is 
K\ = 11. Let us also require the symmetry conditions 

77 77* 77 77* 77 77* 

-^0,1 - -^0,5) ^0,2 - ^o,4> ^0,3 - ^0,3- 

As a result the one-soliton solution simplifies significantly 

Q 10 (x,t) = ^sinh2^ cosh(^(zi + z 2 ))e- i ^ 1 - 22+ ^-^ ) , 

Qn(x,t) = -^^ S inh0 o sinh(^i)e-^ 1+ * 1 \ 



Ai 



IV 



Q 12 (x,t) = ^sinh(20 o )cosh(i/(zi - Z2 )) e - i '*(*i+*a+tfi+fc), 
Q 01 (x,t) = ^^cosh0 o cosh(^ 2 )e-^ 2+ ^), 



where we have used the representation 

7 0,2 = 



F01 = ^sinh^e^ 1 : 



= ^cosh^e^ 2 , 
V2 
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Ai(x, t) = 2 (sinh 2 6q sinh 2 (z^zi) + cosh 2 9q cosh 2 (^£2)) , 



where #0 is an arbitrary real constant. 

If we apply this dressing procedure to the one-soliton solution we obtain a two-soliton solution. 
Iterating this process we can generate multisoliton solutions, i.e. 



Qi s (x,t) -A Q 2s {x,t) 



Qms {.X j i) • 



The particular case s± = s 2 = 1 leads to the compact real form so (5, 0) ~ so (5, R) of the B2- 
algebra. The choice s\ = —s 2 = —1 leads to the noncompact real form so(2, 3) and s± = s 2 = — 1 
gives another noncompact one - the Lorentz algebra so (1,4). If s\ = s 2 = 1 and we identify 

Qoi(x,t) = — Q P oi(x,t), Qw(x,t) = -~E B (x,t), 

% 1 

Qn(x,t) = — E p (x,t), Qn(x,t) = — E a (x,t), 

K K 

then we obtain the system similar to that studied in [TJ [11] which describes Stokes-anti-Stokes 
wave generation. Here Q po \(x,t) is the normalized effective polarization of the medium and 
E p (x, t), E s (x, t) and E a (x, t) are the normalized pump, Stokes and anti-Stokes wave amplitudes 
respectively. In this case the matrix Q(x,t) is defined by 

/ E s E p E a \ 
E* Q pol E a 

E p Qpol Qpol ~ E P 



I 



K Q 



Q{x,t) 

1% 

V E. 

and for 4-wave system we have 



pol 







-e; e* s 



E s 
J 



(Ji - J 2 )E Sit (x,t) - (ii -I 2 )E ByX (x, t) - E p {x,t)Q* pol {x,t) = 0, 
J2Q P oi,t(x,t) - I 2 Q po \,x(x,t) - E*(x,t)E a (x,t) - E p (x,t)E*(x,t) = 0, 
JiE Ptt {x,t) - IiE PtX (x,t) - E a (x,t)Q pol (x,t) + E s (x,t)Q po i(x,t) = 0, 
(Ji + J 2 )E^ t (x, t) - (h + h)E ajX (x, t) + Ep(x, t)Q po i(x, t) = 0. 



Finally we can rewrite it in the form: 

1 0E s (x,t) 0E s (x,t) 



v-i dt 
0Q po i{x,t) 
Ot 



+ 



k-iE p (x, t)Q* ol (x,t) = 0, 



dx 

Kpol 

(E* p (x, t)E a {x, t) + E p (x, t)E*(x, t)) = 0, 



1 0E D (x,t) dEn(x.t) . . ^ . . . . „ 

- g t ' + q x ' ~ K {E & {x,t)Q*p ol {x,t) - E s (x,t)Q pol (x,t)) = 0, 

1 8E a (x,t) 0E a {x,t) 



(3.5) 



v\ Ot 



+ 



Ox 



+ KxE p {x, t)Q po i(x, t) = 0, 



where 



V-l 



dk-i 



h-h 
J1-J2' 



K-i 



h-h' 
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dwi h + h 1 1 

Vl = ITT = 7,7 ' K l = ~ 7 i 7 » K Pol = "I"- 

dfei Ji + J 2 h+h h 

The particular case si = S2 = L J2 = 0, the one-soliton solution obeys (|3.5p . Finally for the 
simplest one soliton solution of (|3.5p the modulus squared of one soliton amplitudes are given by 

2 2 

|£ s O,r,)| 2 = ^-sinh 2 29 cosh 2 z/[( Ji + J 2 )x + Jit], 

g^.2^,2 

| Ji7p (rc, t) 1 2 = — — — sinh 2 #o sinh 2 u( J\x + Jit), 

2,, 2 



|J^a(^,t)| 2 = ^2~ Smh2 2^0 COsh 2 l/[(J! - J 2 )x + Jit], 



1 

g^.2^,2 

|Q P oi(x, t)| 2 = ^ 2 cosh 2 #o cosh 2 vJ 2 x, 

where 

k = — J 2 Ji, Ai(x, t) = 2 (sinh 2 #o sinh 2 (z^( Jix + Jit)) + cosh 2 #o cosh 2 (z/J2x)) . 

Note that the canonical reduction ensures that Ai(x, t) is positive for all x and t, so the solitons 
of this model can have no singularities. 

3.2 Another type of Z 2 reduction 

Let us consider the 7L 2 reduction 

X~(x, t, A) = K 2 { [ X + (x, t, -A)] T } _1 K 2 \ 

where Jf2 = diag (si, s 2 , 1, S4, S5) and si = S5 = ±1, s 2 = S4 = ±1. 
Therefore we have the symmetry conditions 

K 2 U(x,t,-X) T K^ 1 = -U{x,t,X) => Q(x,t) = K 2 Q T (x,t)K^ 1 . 

In particular if we choose Jf2 = 1 then Q(x,t) is a symmetric matrix. 
The invariance condition implies that the dressing matrix gets the form 

( i \\ 1. ■ A ^ K 2 SA{x,t){K2Sr l 
9(x, t, A) = a + 3^ . 

Hence here the poles of the dressing factor form a doublet {A + ,— A + } whose residues are re- 
lated by: 

B(x,t) = -K 2 SA(x,t)(K 2 S)- 1 , A~ = -A+. 
The soliton solution is expressed by the residues of the dressing matrix as follows 

[J,Q(x,t)] = [J,A(x,t) - K 2 SA(x,t)SK 2 }. 
Like before one can present the matrix A(x,t) in the following way 

A(x,t) = X(x,t)F T (x,t), 

where F(x,t) = Sxq(x, t, X + )SFq. In the soliton case it is just the well known plane wave 
F(x,t) = e lX+ ( Jx+It ^Fo. The other factor X(x,t) is given by 

X(x, t) = 2A + (K 2 F - 2\ + SF(F T K 2 F)- l a) (F T K 2 F - A(\+) 2 a(F T XiF)' 1 a) _1 . 
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As a particular case when rankX(x,t) = rankF(x,t) = 1 and a(x,t) = it is simply 

2X+K 2 F(x, t) 
F T (x,t)K 2 F(x,t)' 



X ^ = ^M^TY F-(x,t)K 2 F(x,t) = ^s k e^^Fl k . 



k=l 

Consequently the one-soliton solution reads 

2A+ 



Qij{x,t)=\ FT( x ,t)K 2 F(x,t) VV—^J 'v .o-n— ,o-„, . ^ ,, (3 g) 

0, i = j. 

In the case of B2 algebra there are only 4 independent fields as shown below 

/ Qw{x,t) Q n (x,t) Qi 2 (x,t) 

sis 2 Q w (x,t) Qoi(x,i) Qu(x,t) 

Q(x,t) = sxQn(x,t) s 2 Q i(x,t) Qm{x,t) -Qn(x,t) 

sis 2 Qi 2 (x,t) s 2 Qoi(x,t) Qio(M) 

\ sis 2 Q 12 (x,t) -siQ n (x,t) sis 2 Qio(x,t) 

The corresponding 4-wave system reads 

i(Ji - J 2 )Qio,t(x,t) - i(h - I 2 )Q w ,x(x,t) + ks 2 Q u (x,t)Q 01 (x,t) = 0, 
iJiQu,t(x, t) - iIiQu, x (x, t) + kQ i(x, t)(s 2 Qi 2 (x, t) - Q w {x, t)) = 0, 
i(Ji + J 2 )Qi2,t(x,t) - i(h + h)Qi2,x(x,t) - kQ u (x,t)Qoi(x,t) = 0, 
iJ2Qoi,t(x, t) - H 2 Qoi, x (x, t) + ksiQn(x, t)(Qi 2 (x, t) + s 2 Q 10 (x, t)) = 0. 

Its one-soliton solution when K 2 = 11 is presented by the expressions 

Q10OM) = |^ (e iX+ ^+ J ^+^ +I ^F 0>1 F 0>2 + e - iX+ ^ +J ^ + ^ +I ^F 0A F , 
Qu(M) = (e* + <**+**>fb,i - , 

Q 12 (x,t) = |^ (> + [C*-a>+(^)*]^ 

Q01OM) = (^-+^0,2 + e-K+Vv+'riFoJ , 

where F T (x,t)F(x,t) is obtained from the expression in (|3.6|) putting = 1. After imposing 
the additional restriction Fo t i = Fq^ and Fq i2 = Fq^ we get 

\+ 

Q w (x,t) = sinh 2fl cos A+ [(Ji + J 2 )x + (h + I 2 )t], 

A 2 (x,t) 

2\/2 iA + 

Qn(x,i) = — — r sinh 29q sin A + ( J\x + iii), 

Qi2(ac,*) = ^ sinh 2g cos A+[(Ji - J 2 )x + (h - I 2 )t], 
2\/2A~'~ 

QoiOM) = -1—7 — -- cosh 29q cos X + (J 2 x + I 2 t), 
A 2 (x,t) 



where 



F ,i = sinh 6» , F 0j2 = — j= cosh O , (3.7) 
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#o is a complex parameter and 

A 2 (x,t) = 2(cosh 2 /9 cos 2 X + {J 2 x + I 2 t) - sinh^osin 2 A + ( J x x + ht)). 

In its turn the existence of such reduction leads to the existence of a special class of solitons, 
the so-called breathers, in the case when there are two 7j 2 reductions (canonical one and another 
one of the type mentioned above) applied to the TV-wave systems. 

4 AT-wave equations with a Z 2 X Z 2 -reduction 
and their soliton solutions 

In this section we consider a Z2 x ^-reduced iV-wave system associated to the orthogonal 
algebra. Such systems admit real valued solutions. The soliton solutions can be classified into 
two major types: "doublet" solitons which are related to two imaginary discrete eigenvalues 
(a "doublet") of L and solitons which correspond to four eigenvalues (a "quadruplet"). 

4.1 Doublet solitons 

Let the action of 7L 2 x r L 2 in the space of fundamental solutions of the linear problem is given by 

X -(x,t,X) = K x ((x+^x^A*))" 1 ^ 1 , 
X" t, A) = K 2 (( X + f(x, t, -A)) K~\ 

where K x>2 G SO(n) and [K X ,K 2 ] = 0. Consequently the potential U(x,t,X) satisfies the 
following symmetry conditions 

KxU^x^X*)^ 1 = U(x,t,X), K X TK^ 1 = J, 
K 2 U T {x, t, -X)K^ 1 = -U(x, t, A), K 2 JK 2 - X = J. 

In accordance with what we said in previous chapter the dressing factor g(x,t,X) must be 
invariant under the action of Jj 2 x Z2, i.e. 

K x (g\x, t, X*)\ ^ K^ 1 = g(x,t, A), (4.1) 



K 2 (g T (x,t,-X)) 1 K^ 1 = g(x,t, A). (4.2) 

Let for the sake of simplicity require that K\ = K 2 = 11. As a result we find that the poles of 
the dressing matrix can be purely imaginary, i.e. 

A ± = ±iv, v > 0. 

Thus the invariance condition implies that the dressing matrix gets the form 

g(x,t, X) = H + ^ + SA * ( f[ t)S , A*(x,t) = -A(x,t). 
X — iv X + iv 

Following already discussed procedures we derive in the simplest case that the explicit form of 
A(x,t) is 

A(x,t) = 2 ™ F(x,t)F T (x,t), 

F 1 [x, t)F(x, t) 
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where the vector F(x,t) = e v ( Jx+It >Fo is real. Consequently the soliton solution written in 
a standard matrix notation is 



r 0, 



Qki{x,t) 



if k = L 



2iv 



{F k {x,t)Fi{x,t) 



F T (x,t)F(x,t) 

+ (-l)k+l+i F6 _ k (x,t)F 6 ^(x,t)), ifk^l. 



(4.3) 



In the case of B2 algebra we obtain the following 7L 2 x ^-reduced 4-wave system 

(Ji - J2)qw,t{x,t) - (h - h)qio,x(x,t) + kq u (x,t)q 01 (x,t) = 0, 
J\qii,t{x,t) - hqn,x(x,t) + k(qi 2 {x,t) - q w (x,t))q 01 (x,t) = 0, 
(Jl + J 2 )qi2,t(x,t) - (h + h)qi2,x(x,t) - kq u (x,t)q 01 (x,t) = 0, 
J2qoi,t(x,t) - hqoi,x(x,t) + k(q w (x,t) + qi 2 (x,t))q u (x,t) = 0, 

where qio(x,t), qn(x,t), q\ 2 (x,t) and qoi(x,t) are real valued fields and their indices are asso- 
ciated with the basis of simple roots of B2 introduced in the previous section, i.e. 

Qio(x, t) = iqi (x, t), Qn(x, t) = iqu(x, t), 
Q 12 (x, t) = iq 12 (x, t), Qoi(x, t) = iq i(x, t). 

The "bonding" constant k coincides with the one in the previous examples 

k = J1I2 — J 2 I\- 

Therefore the solution (|4.3p turns into 

, , 2v 
qio[x,t) 

qn{x,t) 

qi2{x,t) 

qoi(x,t) 



F T F 
2v 

Wf 
2v 

Wf 

2v 

Wf 

These solutions can be rewritten in terms of hyperbolic functions as follows 
\v 

q W (x, t) = -^N^ cosh{z4(Ji + J 2 )x + (h + I 3 )t] + 5i + 5 2 }, 



e - v [(Ji+J a )*+(Ji+i a )t] Foili r 0i2 + ^[(^+^+('1+^)^0,5^0,4) , 

e -u[{Ji-J2)x+{h-h)t] p Q x p Q 4 _j_ e v[(Ji-J2)x+(h-h)x] p Q 5 _p Q 2 
e . u{ j 2X+ht) + eU{ j 2X+ht) 



qn(x,t) 



F T F 



iViF , 3 s\nh[u(J lX + ht) + Si], 



qu(x,t) = -jpjT-pNi^ cosh{i/[(Ji - J 2 )x + {h - I 2 )t] + <5i - <5 2 }, 
Av 

ypp-^Foj cosh[^(J 2 2; + I 2 t) + 62], 
F T (x, t)F{x, t) = 2Nl cosh2{v(J lX + ht) +81) + 2N$ cosh2(v(J 2 x + I 2 t) + 5 2 ) +F$ 3 , 



qoi(x,t) 



where we have implied that Fq^ > for k = 1, 2, 4, 5 and therefore the following expressions 

2 -F0,1 

make sense. 



I In A^i = ^0,1*0,5, N 2 = y/F 0t2 Fo A 

z 1*0,2 



16 



V.S. Gerdjikov, N.A. Rostov and T.I. Valchev 



In particular when Fo,i = -Fo,5 and Fq )2 = ^b,4 or in other words S\ = 5 2 = we can apply 
the same representation as in (|3.7|) and as a result we obtain 

q w (x,t) = —— sinh(20 o ) coshi/[(Ji + J 2 )x + (h + I 2 )t], 

qn(x,t) = smh.6osmhu(Jix + lit), 

qi2(x,t) = —— sinh(26>o) coshz/[(Ji - J 2 )x + {I\ - I 2 )t], 

qoi{x,t) = — cosh 9q cosh v{J 2 x + I 2 t), 

where 9o and 

Ad(x, t) =2 (sinh 2 9q sinh 2 v{J\x + I\t) + cosh 2 6q cosh 2 u(J 2 x + I 2 t)) . 
4.2 Quadruplet solitons 

There is another way to ensure the Z* 2 x Z 2 invariance of the dressing factor. This time we 
consider dressing factors g(x, t, A) with two more poles. The requirements (|4.1|) - (|4.2|) lead to 
the following dressing matrix 

A(x,t) , K^A^x^i^S)- 1 KzSAfatyKzS)- 1 

g(x,t, A) — 1L 



A-A+ A-(A+)* A + A- 

A + (A+)* 



(4.4) 



By taking the limit A — > 00 in equation (|2.5p and taking into account the explicit formula (|4.4p 
one can derive in the soliton case Qo(x,t) = the following relation 

[J,Q](x,t) = [J,A + K 1 SA*SK 1 -K 2 SASK 2 -K 1 K 2 A*K 2 K 1 ](x,t). (4.5) 

Like in previous considerations we decompose the matrix A(x, t) using two matrix factors X(x, t) 
and F{x, t) and derive some differential equation for F(x, t) which leads to 

F(x,t) = e iX+ ^ Jx+I ^F . 

The linear system for X(x,t) in this case is following 

/ K 2 SA(x,t)SK 2 KjSA^x^SKi K&A^xrfKzKA 

V 2A+ 2iv 2(i J [Xl ' 

= X(x,t)a(x,t), (4.6) 

where a(x,t) is a linear function of x and t as follows 

t) = iF^(Jx + It)SF . 



Starting from the equation (14. 6j) by multiplying with K\S, K 2 S and performing complex con- 
jugation if necessary, we can derive the following auxiliary linear system 

SF(x, t)=Xa + - Z^— + W-^- , 

SG(x, t) = X— + Ya + Z— - W^—, 
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SH(x , t) = X <£M1 + r <^> + Za « + 



2iu 



2/' 



2(AH 



2/x 2if 
where we introduced auxiliary entities 



2(A+)* 



Y(x,t) = K 2 SX(x,t), Z(x,t) = K 1 SX*{x,t), W(x, t) = K 1 K 2 X*(x, t), 
G(x,t) = K 2 SF(x,t), H(x,t) = K-lSF^x^), N(x,t) = K 1 K 2 F*{x,t), 
(F, H) = F T SH. 



In matrix notations this system reads 



( SF, SG, SH, SN ) = ( X, Y, Z, W ) 



( a a b c \ 

a a c b 

b* c* a* a* 

\ c* b* a* a* ) 



where 



_ (F{x,t),G{x,t)) (F(x,t),H(x,t)) 

Q(X,l) — „, 1 ) 0{X,l) — 



2A+ 



2iu 



c(x, t) 



(F(x,t),N(x,t)) 
2/i 



To calculate X(x,t) we just have to find the inverse matrix of the block matrix shown above. 
In the simplest case when rankX(x,t) = rankF(x,t) = 1 and a = we have 



/ x \ 




/ a* 


b 


—c 


Y 


1 


a* 


—c 


b 


Z 


A(x,t) 


-b -c 





a 


V w J 




\ -c -b 


a 






\ 



( SF \ 

SG 
SH 
V SN ) 



where 

A(x, t) = \a(x, t)\ 2 + b 2 (x, t) - c 2 . 
Finally putting the result for X(x,t) in (|4.5jl we obtain the quadruplet solution 

Q(x,t) = — [(a*(x,t)K 2 F(x,t) +b(x,t)K 1 F*(x,t) - c(x,t)K 1 K 2 SF* (x,t)) F T (x,t) 

- K 2 S (a*(x, t)K 2 F(x, t) + b(x, t)K 1 F*{x, t) - c(x, t)K 1 K 2 SF*(x, t)) F T (x, t)SK 2 
+ KiS (a(x, t)K 2 F*(x, t) - b(x, t)KiF{x, t) - c(x, t)KiK 2 SF{x, t)) F ] (x, t)SK\ 

- K X K 2 {a{x, t)K 2 F*(x, t) - b(x, t)K x F(x, t) - c(x, t)K 1 K 2 SF(x, t)) F f (x, t)K 2 K x . 

Consider the 4-wave system associated with the B 2 algebra. Let K\ = K 2 = 1 then its 
generic quadruplet (or breather-like) solution is 



9io(x,i) = — Im 
4 

qn(x,t) = — Im 
4 

qi 2 {x,t) = — Im 



a*Ni cosh(y?i + ip 2 ) 



imN 



[IV 



- (//cosh(<£* + ip 2 ) - iucosh(ipl - ip 2 )) 



a*Ni sinh(^i) 



■iV*sinh(^) 



Ft 



0,3) 



a*Ni cosh(( / 9i — (p 2 ) - (/icosh(</?* — ip 2 ) — iv cosh(y?* + ip 2 )) 



N 2 , 
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q i(x,t) = — Im 



where 



a(x, t) 
b(x,t) 



1 



H + iv 

m(x, t) 
iv 



a*iV 2 cosh(y? 2 ) Nn coshfwol 



F 2 

iVf cosh 2v?i + iV| cosh 2<p 2 + 



c(x, t) 



777,(3;, i) 



m(x, t) = |Aq| 2 cosh(2Re <pi) + |iV 2 | 2 cosh(2Re <p 2 ) + 



1 . F ( 



(fi a (x, t) = i\ + (J a X + I a t) + - log 



O.o- 



F, 



0,3) 



1-^0,3 1 2 

2 ' 



2 ^0,6-0- 



a = 1,2. 



^cr — V -Pb,o--f0,6-cr, 



5 Conclusion 

New N-w&ve type equations related to the orthogonal algebras and obtained from the generic 
ones via additional reductions are analyzed. These new nonlinear equations are solvable by the 
inverse scattering method. In particular, the systems related to the 50 (5) algebra involve 4 waves. 
Imposing Z 2 and Z 2 x Z 2 reductions on them we obtain 4-wave systems of physical importance. 
The 4-wave system with a canonical Z 2 reduction describes Stokes-anti-Stokes Raman scattering. 
The Z 2 x Z 2 -reduced N wave system possesses real valued solutions while other types of N wave 
equations have complex valued solutions. This determines its mathematical importance. 

The soliton solutions of these integrable systems are parametrized by two types of parameters: 
by the discrete eigenvalues of L and by the "polarization" vector Fq. Therefore the problem 
of classifying all types of one-soliton solutions of the N-w&ve equations is equivalent to that 
of classifying all possible types of discrete eigenvalues and polarization vectors. When both 
reductions are applied we have two different configurations of eigenvalues for the Lax operator L: 
doublets and quadruplets. This situation is analogous to one encountered in the sine-Gordon 
case, which was the reason to call quadruplet solitons 'breather'-like solutions |12j . 
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